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Abstract 

We give exact and approximation algorithms for computing the Gromov hyperboHcity of an 
n-point discrete metric space. We observe that computing the Gromov hyperbolicity from a 
fixed base-point reduces to a (max,min) matrix product. Hence, using the (max,min) matrix 
product algorithm by Duan and Pettie, the fixed base-point hyperbolicity can be determined 
in 0(n^'^®) time. It follows that the Gromov hyperbolicity can be computed in 0{n^^^^) time, 
and a 2-approximation can be found in 0{n^'^^) time. We also give a (2 logj n)-approximation 
algorithm that runs in OirP') time, based on a tree-metric embedding by Gromov. 

1 Introduction 

Gromov introduced a notion of metric-space hyperbolicity [H [6] using a simple four point condition. 
(See Section fl.ll ) This definition is very attractive from a computer scientist point of view as the 
hyperbolicity of a finite metric space can be easily computed by brute force, by simply checking 
the four point condition at each quadruple of points. However, this approach takes 0(n^) time for 
an n-point metric space, which makes it impractical for some applications to networking. Knowing 
the hyperbolicity is important, as the running time and space requirements of previous algorithms 
designed for Gromov hyperbolic spaces are often analyzed in terms of their Gromov hyperbolicity [31 
[7]. So far, it seems that no better algorithm than brute force was known for computing the Gromov 
hyperbolicity [2]. In this note, we give faster exact and approximation algorithms based on previous 
work on (max-min) matrix products by Duan and Pettie 0], and the tree- metric embedding by 
Gromov |6j. 

The exponent of matrix multiplication ^ is the infimum of the real numbers a; > such that two 
n X n real matrices can be multiplied in 0{n^) time, exact arithmetic operations being performed 
in one step [8]. Currently, ^ is known to be less than 2.3727 [9]. In the following, a; is a real number 
such that we can multiply two n x n real matrices in 0{n^) time. 

Our algorithm for computing the Gromov hyperbolicity runs in 0(n(^+'^)/^) time, which is 
0{n^'^^). (See Section [2.11 ) For a fixed base-point, this improves to 0{n^^^^^/'^), which also yields 
a 2-factor approximation for the general case within the same time bound. (See Section 12.21 ) 
Finally, we give a quadratic-time (2 log2 n)-approximation algorithm. (See Section [2?3l ) 
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1.1 Gromov hyperbolic spaces 

An introduction to Gromov hyperbolic spaces can be found in the article by Bonk and Schramm [1], 
and in the book by Ghys and de la Harpe [5]. Here we briefly present some definitions and facts 
that will be needed in this note. 

A metric space (M, d) is said to be (5-hyperbolic for some (5 ^ if it obeys the so-called four 
point condition: For any x,y,z,t G M, the largest two distance sums among d{x,y) + d{z,t), 
d{x,z) + d{y,t), and d{x,t) + d{y,z), differ by at most 26. The Gromov hyperbolicity 5* of {M,d) 
is the smallest 5* ^ such that (M, d) is (5*-hyperbolic. 

For any x,y,r G M, the Gromov product of x,y at r is defined as 

ix\y)r = ^ {d{x, r) + d{y, r) - d{x, y)) . 

The point r is called the base point. Gromov hyperbolicity can also be defined in terms of the 
Gromov product, instead of the four point condition above. The two definitions are equivalent, 
with same values of 5 and 5*. So a metric space {M,d) is 5-hyperbolic if and only if, for any 
x,y, z,r £ M 

{x\z)r ^ min{(x|y),., {y\z)r} - 6. 

The Gromov hyperbolicity 6* is the smallest value of 6 that satisfies the above property. In other 
words, 

5* = max {mm{{x\y)r, {y\z)r} - {x\z)r} . 

x,y,z,r 

The hyperbolicity 5r at base point r is defined as 

6r = max {mm{{x\y)r, iy\z)r} - {x\z)r} . (1) 

x,y,z 

Hence, we have 

6* = max5r. (2) 

r 

2 Algorithms 

In this section, we consider a discrete metric space (M, d) with n elements, that we denote xi, . . . ,Xn. 
Our goal is to compute exactly, or approximately, its hyperbolicity 6* , or its hyperbolicity 6r at a 
base point r. 

2.1 Exact algorithms 

The (max,min)-product A0 B two real matrices A, B is defined as follows: 

{A (g) B)ij = maxmin{Ajfc, Bkj}. 

k 

Duan and Pettie gave an 0(n*-'^^'^''/^)-time algorithm for computing the (max,min)-product 
of two n X n matrices. 

Let r be a fixed base-point. By Equation ([T]), if ^ is the matrix defined by Aij = {xi\xj)r for any 
i,j, then 6r is simply the largest coefficient in {A^ A) — A. So we can compute 6r in 0{n^^~^^^^'^) 
time. Maximizing over all values of r, we can compute the hyperbolicity 6* in 0{n^^^^^^'^) time, by 
Equation ([2|). 
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2.2 Factor-2 approximation 

The hyperbolicity 5r with respect to any base-point is known to be a 2-approximation of the 
hyperbohcity 6* [1]. More precisely, we have dr ^ 6* ^ 26r- So, using the algorithm of Section [2?T| 
we can pick an arbitrary base-point r and compute 6r in 0(n('^+'^)/-^) time, which gives us a 2- 
approximation of 6*. 

2.3 0(logn)-factor approximation 

Gromov [6] (see also the article by Chepoi et al. [3l Theorem 1] and the book by Ghys and de la 
Harpe [5) Chapter 2]) showed that any 5-hyperbolic metric space {M,d) can be embedded into a 
weighted tree T with an additive error 2(51og2ra, and this tree can be constructed in time 0{'n?). 
In particular, if we denote by dx the metric corresponding to such a tree T, then 

d(a, h) — 25* logg n ^ dT{a^ b) ^ d{a, h) for any a,b ^ M. (3) 

This construction can be performed without prior knowledge of 5*. 

We compute D = max^ fegjv/ d{a, b) — dxia, b) in time O(n^). We claim that: 

6* i^D 26* log2 n. (4) 

So we obtain a (2 log2 ?T-)-approximation 2D of 5* in time 0(n-^). 

We still need to prove the double inequality It follows from Equation ([3]) that d(a,b) — 
dT{a,b) ^ 2(^*log2ra for any a,b, and thus D ^ 25*log2n. In the following, we prove the other 
inequality. 

For any x, y, z, t, we denote by 5{x, y, z, t) the difference between the two largest distance sums 
among y) + d{z, t), d{x, z) + d{y, t), and d{x, t) + d{y, z). Thus, if for instance d{x^ y) + d{z, t) ^ 
d{x, z) + d{y, t) ^ d{x, t) + d{y, z), we have S{x, y, z, t) = d(x, y) + d(z, t) — d{x, z) — d{y, t). We also 
need to introduce the difference 5t{x, y, z, t) between the two largest sums among drix, y) + dT{z, t), 
drix.z) +dT{y,t), and dT{x,t) + dT{y,z). 

For any a,b G M, we have d{a, b) — D ^ dxia, b) ^ d{a, b), so 6{x, y, z, t) — 6t{x, y, z, t) ^ 2D, 
because in the worst case, the largest sum with respect to d is the same as the largest sum with 
respect to dx, and the second largest sum with respect to dr is equal to the second largest sum 
with respect to d minus 2D. But by construction, dx is a tree metric [3], so 5T{x,y,z,t) = for 
any x, y, z, t. Therefore 6{x, y, z, t) ^ 2D for any x, y, z, t, which means that 6* ^ D. 
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